Array- induced collective transport in the Brownian motion of 
coupled nonlinear oscillator systems 



Zhigang Zheng 1,3 Bambi Hu 1 ' 4 and Gang Hu 2,3 
1 Department of Physics and Center for Nonlinear Studies, Hong Kong Baptist University, 

Kowloon Tong, Hong Kong, China 
2 Center of Theoretical Physics, Chinese Center of Advanced Science and Technology (World 

Laboratory), Beijing 8730, China 
3 Department of Physics, Beijing Normal University, Beijing 100875, China 
4 Department of Physics, University of Houston, Houston TX 77204, USA 



Abstract 

Brownian motion of an array of harmonically coupled particles subject to 
a periodic substrate potential and driven by an external bias is investigated. 
In the linear response limit (small bias), the coupling between particles may 
enhance the diffusion process, depending on the competition between the 
harmonic chain and the substrate potential. An analytical formula of the 
diffusion rate for the single-particle case is also obtained. In the nonlinear 
response regime, the moving kink may become phase-locked to its radiated 
phonon waves, hence the mobility of the chain may decrease as one increases 
the external force. 
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I. INTRODUCTION 



Brownian motion of particles subject to periodic substrate potentials and external forces 
gained great interest due to its wide applications and practical importance in connection with 
transport processes in many fields including damped pendula, superionic conductor, Joseph- 
son tunneling junction, vortex motion in high-T c oxide superconductors, phase-locked loop, 
rotation of dipoles, charge- density wave, dislocation, and so on [1-7]. Although the collective 
transport of coupled nonlinear oscillators has been studied recently for the two-dimensional 
case in relation to studies on adsorbate islands and monolayer films on surfaces [8,15,17], 
explorations of collective behaviors for the one-dimensional coupled nonlinear oscillators are 
still fundamentally important in understanding many physical systems. In dimensionless 
form, the Langevin equation in describing the ID case might be written as 



x 3 +7 X] +dsin (^p) = J2 + m + (1) 



where Xj represents the coordinate of the i-th particle and Xj= dxj/dt its corresponding 
velocity, 7 is the friction coefficient, d is the height of the periodic potential, and b denotes 
the period of the substrate potential. V(xj,Xi) reflects the interaction between the j-th and 
the i-th particle. Fj(t) and £j(t) denote the external driving force and the thermal fluctuation 
induced random force on the j-th particle, respectively. In this paper we focus on the case 
when the external driving is uniformly constant, i.e., Fj(t) = F. The thermal noise is 
frequently assumed to be a both spatially and temporally uncorrelated Gaussian-type one: 

&(*)> = 0, <&(*)fc(0> = 2 7 *bTM(* - 0, (2) 

where ks is the Boltzmann constant, and T is the environmental temperature. When the 
interaction among particles is very weak, in many cases one may approximately treat the 
above problem in terms of Brownian motion of a single particle under a biased periodic 
potential. The corresponding Fokker-Planck equation then reads [1] 

dW . dW d ( . Tir , rT1 dW\ 

Here W(x, x, t) denotes the probability distribution of the particle, where the indices of par- 
ticles are omitted. This equation is sufficiently complicated to be analytically solved so that 
a closed analytical solution in describing the diffusion process is far from obtainable. Based 
on the inverse-friction expansion method, we recently gave an analytical perturbative solu- 
tion of the mobility for the single-particle case valid for arbitrary friction cases [9]. When the 
interaction among particles can no longer be ignored, one has to treat the coupled case (1). 
Interactions between particles introduce new time scales, which leads to more complicated 
phenomena, and furthermore an analytical treatment of the corresponding Fokker-Planck 
equation is almost impossible [3]. One of the simplest models in describing the competition 
between the coupling and the substrate potential is the well-known Frenkel-Kontorova sys- 
tem [10], which describes a chain of particles with nearest-neighboring harmonic couplings 
and subject to a periodic potential, where the interaction is reduced to 
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V(xj, Xj-x) = l -K[xj - Xj-! - a) 2 , (4) 

where K and a are the coupling strength and the static length of the spring, respectively. 
Now there are two competing lengths: a and b. The winding number (or the frustration) 
is defined as 5 = b/a, which may strongly affect the spatial configuration of the system. 
During the past few years, the FK model was applied to investigations of the ground state 
of competing systems, and the commensurate-incommensurate (CI) phase transitions were 
found and theoretically explored [11]. The theory developed by Aubry stands as one of 
the deepest achievements in theoretical comprehension of the physics of modulated phases 
[12]. Dynamics of the FK chain was also explored in relating to many fields, such as charge- 
density wave (CDW), nano-tribology and surface problem, self-organized criticality (SOC) 
and Josephson-j unction arrays (JJA) and ladders (JJL) [13, 14]. In the following investiga- 
tions, we simply set b to be 2n. Then the Langevin equation connected to the FK chain can 
be followed from (1): 

Xj +7 Xj +dsinxj = K(xj + i — 2xj + a;j-i) + F + £j(t). (5) 

One may notice that the static length of the spring a does not enter (5), but it may play 
a very significant role in describing the motion of the coupled systems. The corresponding 
Fokker-Planck equation reads 



dW _ v 
dt ~^ 
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dW dU dW d (. ,„ , 0W" 



x, — — + — — + 77 — \x j W + k B T 



1 dxj dxj dxj d Xj \ 3 ' '" D ~ dx 3j j 



(6) 



where W = W({xj} , {ij} ,t) is the joint probability distribution function, and here we use 
the total potential U({xj}), and for the case of the FK model U({xj}) = Y^[d(l — cosxj) + 
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^K(xj + i — Xj — a) 2 }. In the high-friction limit one may obtain the Smoluchowski equation 
by averaging velocities, which can be written as 

dt 7 ? dxj i^dxj ^ B dx^j ' ^ 
where P({xj} ,t) = J W({xj} , {ij} , t) TT d Xj is the reduced probability distribution for only 
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spatial variables. In many cases one is interested in the mobility of the chain: 

(8) 

where < . > includes both time and particle averages. This definition was introduced in 
discussions of Brownian motion of a single particle in the biased periodic potential. We 
will give a unified solution for the single particle case in the following discussion. Ana- 
lytical investigations of the coupled case may be more difficult, therefore we will discuss 
this problem mainly in terms of numerical simulations. The collective diffusion coefficient, 
which is relevant for studying commensurability effects [15], is given by the linear part of 
the mean-square displacement: 
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(9) 



This coefficient is connected to the mobility in the small force limit (linear response) [1] by 

D(T,K,8) = k B T]imn(F,T,K,8). (10) 

F — »0 

One may notice that the diffusion coefficient relates to the temperature, the coupling strength 
and the frustration. In particular, the dependence on the coupling strength and the frus- 
tration does not happen for the single particle case. Studies of their dependences are very 
interesting and significant. In the following discussions, we will investigate both the linear 
and the nonlinear response regimes. We find that in the limit of linear response, the coupling 
between particles can greatly increase the diffusion rate, which may be important in real- 
istic experiments and technologies. Additionally, we shall investigate the commensurability 
effects. In the vicinity of the golden mean 8 = 8q = (y/5 — l)/2, the chain possesses the 
maximum diffusion coefficient. In the nonlinear response regime, the transport is dominated 
by the strong resonant behavior due to the competitive phase-locking between the travel- 
ling wave and its radiated linear phonons. This resonance may lead to the suppression of 
the mobility. In the following investigations, we mainly perform numerical simulations for 
the coupled case. The fourth-order Runge-Kutta integration algorithm is applied and the 
time step is adjusted according to the numerical accuracy. Periodic boundary conditions are 
applied, i.e., Xj + ^(t) = Xj(t) + 2tvM, where M is an integer that counts the net number of 
kinks trapped in the ring, therefore the frustration is 8 — M/N and the spring constant will 
be a = 2ir5. Throughout the paper 7 = 0.1 and d is set to be 1. 



II. LINEAR RESPONSE: ENHANCEMENT OF THE DIFFUSION 

In Fig. 1(a), we give the relation between the diffusion coefficient D and the coupling 
strength K for the incommensurate case (e.g., the gold mean case 8 = 8 a) and for different 
temperatures. The first phenomenon we observe is that the relation between the diffusion 
coefficient and the coupling strength is not monotonic. For weak coupling strengths, the 
diffusion process is suppressed; When K > K cl (the first threshold), the diffusion coefficient 
begins to increase as one increases the coupling strength. Near this critical value, the 
diffusion coefficient D is found to obey the following relation: 

D = D (T)\K-K cl \ a(T) , (11) 

where the scaling exponent a(T) decreases with increasing the temperature T. At a second 
threshold K c2 , D begins exceeding the single-particle value (K — 0, the non- interacting 
case), exhibiting an array- enhanced diffusion process. For the moderate coupling strength, 
we find that the D — K relation obeys the power law: 

D oc K^ T \ (12) 

where the scaling exponent scales with the temperature as (3(T) ~ T -1 / 2 , i.e., it decreases 
with increasing the temperature. At the high coupling constant, the diffusions coefficient 
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may saturate to a value much higher than the single-particle case. This behavior is very 
interesting, because one must overcome a threshold coupling before one can get a higher 
diffusion rate. Additionally K cX and K c2 increases with the temperature, i.e., for higher 
temperatures, one must overcome stronger coupling thresholds to get higher diffusion rates. 
The result indicates that if one introduces some coupling between the particles (elements) 
then a higher diffusion rate than uncoupled systems can be achieved. In many realistic 
applications one hopes that the diffusion process can be improved as fast as possible. Our 
exploration indicates that for the incommensurate case the coupling between particles may 
enhance the diffusion process. 

The above behavior can be heuristically interpreted, which is a typical consequence of the 
competition between order and disorder. The mechanism of order comes from the coupling 
among particles, where the coupling tends to organize the chain to move in a collective 
way; the mechanism of disorder is the thermal fluctuation, which tends to destruct the 
ordered motion. For very weak couplings, the thermal noise may dominate, thus particles 
cannot organize themselves very well to diffuse collectively. In this case interactions between 
elements introduce another resource of the dissipation that leads to higher friction [18]. 
When the coupling between particles increases and exceeds a threshold, then the role of 
disorder may be supressed, and particles can gradually move collectively, leading to high 
diffusion rates. At higher temperatures, the noise produces much more disorder, hence the 
chain needs a stronger coupling to organize the collective diffusion. 

All the above discussions are valid for the incommensurate case. In Fig. 1(b), we give the 
numerical result of the diffusion coefficient against the winding number (frustration) 8 for 
different temperatures, and the coupling strength K = 1.0, which strongly enough for the 
chain to diffuse. It can be found that for high temperatures, the commensurability effect is 
not very significant. Because the curve is symmetric about 1/2, we only discuss the range 
5 G [0, 1/2]. For small winding numbers, the relation is a linear one (see Fig. 1(b)): 

D = D T5, (13) 

where the slope is proportional to the temperature. This indicates that for small frustrations, 
the diffusion is suppressed by introducing the coupling among particles. The curve becomes 
flat when 5 further increases. This result indicates that although the particles are coupled, 
for small frustrations, the diffusion process is still very slow. In this case couplings between 
particles act as an additional source of dissipation [18]. This suppression of the diffusion is 
because the PN barrier is high (for the coupled case the motion of the chain is dominated 
by the moving kink, as pointed out below. The Peierls-Nabarro barrier corresponds to the 
barrier for the kink translation along the chain. It is also the minimum energy necessary to 
move the kink along the chain.), and the collision of particles with the substrate potential 
plays a more significant role, thus one needs a higher activated energy to overcome the PN 
barrier. Commensurate effect is clearly shown when the temperature decreases. For the case 
ksT = 0.25, we observe that the diffusion rate reaches a maximum value at approximately 
5 = 1/3 and 2/3. These two values are very close to the golden mean value 8q = (V§ — l)/2 
and 1 — 5c- In order to get a higher diffusion, one should choose winding numbers near 
the golden mean. This result can be applied to the dynamics of Josephson junction arrays, 
where the frustration can be altered by changing the magnetic field strength [19]. 
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The relation between D and T is not trivial. We first discuss the non-interacting case. For 
the very small damping constant and high enough temperature case, the diffusion coefficient 
was approximately obtained as [1] 

„ nk B T . 2d . 

D=— — exp(-- , 14 

2 7 PV k B T h v ; 

i.e., the diffusion rate increases with the temperature. In the high damping and low tem- 
perature limit, the diffusion is dominated by thermally activated hoppings, then one has the 
following formula [1]: 



D = — 

7 



I< d 



k B T' 



-. 9 

2nd . 2d ^ x 

^— exp( -^»- (15) 



where I n (x) is the modified Bessel function, one may find the same Boltzmann factor (Ar- 
rhenius form) exp(— 2d/k B T) as the low-damping case, except the prefactor. The common 
factor comes from the thermal-fluctuation- induced hopping effect. In general cases, we are 
still able to derive a unified formula of the diffusion rate. We investigated the mobility and 
got a successively perturbative solution of the mobility [9]. By using the relation between 
the mobility and the diffusion coefficient (10), in the linear response limit, F — > 0, we have 
the formula: 

(2ir) 2 k B T 

~ 7 [n(d,0)A(d,0)-A(7,T,d,0)]' (16) 
where Q(d, F), A(d, F) and A(7, T, d, F) are given by 

Q(d, F) = r exp[/(x) - Fx}dx, 
Jo 

A(d, F) = / exp[-f(x) + Fx]dx, 
Jo 

A( 7 ,T,rf,F) = f 2n dxex P [-f(x) + Fx] f q(0 exp[/(0 - F£]dt, (17) 
jo Jo 



where f(x) = —dcosx. The kernel function q(x) can be obtained from the following operation 
(see [9] for a detailed discussion): 

q(x) = ^7 H c, (18) 

A 

where c is a constant and H is a continued-fraction operator acting on the constant c: 



A A 1 A + 1 A 

H=/ -^^ a ; A+ ~ K » (19) 

I-^sK ■ ■ ■ K 



A A + A ~ 

where I is the unit operator and K , K are given as 



A + r, — ^ 9 
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A - A + 

K —K 4 



1 \df(x) 



-F 



(20) 



^/k^T I dx 



This is a unified formula of the diffusion coefficient valid for arbitrary damping and temper- 
ature cases. 

When the particles are coupled in an array way, they will move collectively, leading to 
the wave motion. In this case, the diffusion is mainly related to the kink motion (travelling 
wave) [14, 16]. The kink describes the minimally possible, topologically stable compression 
of the commensurate structure. The kink is a quasiparticle, characterized by an effective 
mass, rest energy and the Peierls-Nabarro (PN) amplitude E PN . In the low temperature 
case, the diffusion rate can be given as [17] 



where D is the prefactor that scales with E PN as D oc E P ' N for the low-damping case and 
D oc E PN for the high-damping case. In the strong coupling case, one may have 



Eq. (22) indicates that there are two competing terms. For very high temperatures, the 
diffusion rate may be dominated by D . In the low-temperature regime, the two terms may 
compete and the exponential term plays an important role. Because the PN barrier closely 
relates to the commensurability effect, the diffusion rate will be strongly affected by the 
commensurability of the system. It may be easily found that for both the single particle and 
the coupled particle systems, an Arrhenius factor always exists, except that for the single 
case the factor relates to the potential barrier, but for the chain it relates to the PN barrier. 

III. NONLINEAR RESPONSE: SUPPRESSION OF THE MOBILITY 

When one increases the external driving force, the transport behavior becomes a drift 
one. In this case one usually uses the mobility ji = (v) / F to describe the transport process 
of the system. In Fig. 2(a) and (b), we give the the numerical /i — F relation for 5 = 3/8 
and 1/8, respectively. The mobility of the single particle for k B T = 0.50 is also plotted to 
make a clearer comparison. The first result one may clearly observe is that when the driving 
force is small, the mobility of the chain is much higher than that of the single-particle case. 
As one increases the force, the mobility decreases. These two phenomena are much different 
from the single-particle case (see the single-particle line in Fig. 2). Another novel behavior 
different from the single-particle case is that at the left side of the decreasing lines(e.g., 
different lines between 0.15 and 0.4 in Fig. 2(a) corresponding to different temperatures, 
see the left side F = 0.15), the mobility decreases when one increases the temperature 
(see, for example, F = 0.15 in the inset of Fig. 2 (a), the mobilities corresponding to higher 
temperatures are smaller.). In the middle of the decreasing lines (e.g., F = 0.25 in Fig. 2(a)), 
the mobility is almost unaffected by the increasing temperature. At the right hand of the 



D = D exp( 




(21) 




(22) 



7 



decreasing lines, the mobility increases with increasing the temperature, which is a natural 
consequence (Similar phenomenon could also occur at approximately F = 0.1). The above 
anomalous phenomena are consequences of the interaction among particles. In fact, one may 
study the noiseless case to find this intrinsic reason. In this case it is found [20-22] that the 
motion of the chain is dominated by the moving localized kink, i.e., a distorted travelling 
wave that is composed of a moving kink and the oscillating linear wave around it. In the 
underdamped case, the moving kink may become phase-locked to its radiated phonons [14]. 
We have obtained a mean-field formula of the resonance-velocity spectrum (see [21] for a 
detailed derivation): 



v( mi ,m 2 ) = ^Jp + 4K sm 2 (^^). (23) 
mi V rrii 

The resonance is denoted by a pair of integers (mi,m 2 ). Here (3 =< J2f=i jf cos(x*(t)) > is 
called the contraction factor, which describes the collective effect of the chain (kink). Here 
x*(t) denotes the steady state of the system (5), and < . > represents the time average. 
In Fig. 3 (a), we give the relation between (3 and 5 for F — 0. In this case, commensurate 
effect can be clearly observed, where some rational peaks correspond to commensurate cases 
5 = 0,1,1/2,1/3,2/3,.... In Fig. 3(b), we plot the contraction factor (3 vs. the external 
driving F . The dashed horizontal line is the value for F = 0. It is clearly illustrated that (3 
varies with F in a non-monotonic way. The value of j3 for F ^ is higher than that of F = 0, 
i.e., the chain is further contracted when one increases F. At F pa 0.47, (3 suddenly drops 
off to a negative value. The compressed topology of the chain (kink) becomes extended, 
indicating a transition to a new state. This transition is also shown in Fig. 2(a) and (b), 
where the mobility jumps to a much higher value. We called this regime with a high mobility 
the "high- velocity regime". This effect will be described in another work, which is beyond 
our present interest. 

The resonance behavior between the travelling wave and its radiated phonons remains 
robust even when finite temperature effect (thermal noise) is considered. The mode locking 
behavior leads to resonant steps of the averaged velocity of the chain when one adiabatically 
increases the external force, thus a decrease of the mobility with increasing forces can be 
expected. When the driving force is increased along a resonant step, the input energy is 
mainly consumed for amplifying the excited linear phonon and even exciting new linear 
modes, while the average velocity of the chain remains nearly constant. Thus the coupling 
between particles here acts as an additional source of dissipation. Due to the amplification of 
the linear wave, the resonance will be gradually unstable and eventually destructed when one 
increases the external force. As one increases the temperature, the excited linear waves will 
be further amplified. At the right end of a resonant step these very strong linear waves would 
destruct the resonances, thus the mobility increases. This is the so-called " smearing effecf . 
At the left right hand of the resonance, the noise can only slightly affect the resonance, hence 
the coupling between particles may stabilize the resonance. In this case the noise tends to 
drive the resonance state to a lower resonance state, thus the mobility decreases when the 
temperature is increased. This interprets the observation reported above. To get a deeper 
understanding, in Fig.4, we give the relation between the mobility and the temperature for 
different external forces F = 0.05 and 0.10 at 5 = 3/8. F = 0.10 corresponds approximately 
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to the left hand of the resonance and F = 0.05 the middle of the resonance. It can be clearly 
shown that in the middle of the resonance, the mobility keeps nearly constant; while at the 
left hand, the mobility decreases with the temperature. 

In conclusion, in this paper we investigated the collective transport behavior of the 
underdamped Frenkel-Kontorova chain under a constant external driving force and the in- 
fluence of the environmental noise, and found complicated behaviors. Results reported in 
this paper should be valuable of applications in many physical cases, such as charge-density 
waves, Josephson junction arrays and ladders and tribology. Some results need a theoretical 
exploration, and studies along this line are currently in progress. 
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FIGURES 
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FIG. 1. (a). The diffusion coefficient D of the FK chain versus the coupling strength K 
for different temperatures fegT = 0.25, 0.5 and 1.0. The vertical coordinate is shown by using a 
logarithm scale. The diffusion coefficient first decreases against the coupling and then increases 
and exceeds the value for the uncoupled case, indicating a competition between order and disorder. 
For moderate coupling, the D — K relation is a power type. (b). The diffusion rate D against the 
winding number (frustration) 5 of the chain for different temperatures fc^T = 0.25, 0.5 and 1.0. 
The diffusion indicates a strong commensur ability effect. 
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FIG. 2. The mobility of the chain \i varying with the external driving force F for cases (a): 
N = 8, M = 3 and (b): N = 8, M = 1 for different temperatures. The comparison line corresponds 
to the case of a single particle at the temperature fcgT = 0.5. The inset of (a) is plotted to make 
a clearer observation. In many regions the mobility decreases with increasing the external force 
for the coupled case. This anomalous behavior is the consequence of the competition between 
the moving kink (traveling wave) and the phonons radiated by the collision between the harmonic 
chain and the periodic potential. 
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FIG. 3. (a): The contraction factor (3 varies with the frustration 5 with F = for K=0.5 
and 1.0. Peaks at the rational 5 = 0, 1, 1/2, ... shows the commensurate effect, (b): (3 against the 
external force F with 8 = 3/8. The dashed horizontal line corresponds to the value of (3 at F = 0. 
The hysteresis loop can be observed due to the bistability of the system. 




FIG. 4. The mobility of the chain fj, varying with the the temperature ksT for the case N = 8, 
M = 3 for different external drives F = 0.05 and 0.1. The mobility remains almost unchanged or 
decreases, which is much different from the single-particle case. 
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